Collective Pinning of a Frozen Vortex Liquid in Ultrathin Superconducting 

YBaaCusOr Films 
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The linear dynamic response of the two-dimensional (2D) vortex medium in ultrathin YBa2Cu307 
films was studied by measuring their ac sheet impedance Z over a broad range of frequencies uu. 
With decreasing temperature the dissipative component of Z exhibits, at a temperature T*{uj) well 
above the melting temperature of a 2D vortex crystal, a crossover from a thermally activated regime 
involving single vortices to a regime where the response has features consistent with a description 
in terms of a collectively pinned vortex manifold. This suggests the idea of a vortex liquid which, 
below T*{uj), appears to be frozen at the time scales 1/uj of the experiments. 
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It is generally accepted that the effect of weak disorder 
on flux lattices in superconductors results in some sort 
of glassy state MM- In a rigorous sense, however, the 
existence of a truly superconducting vortex glass phase, 
such that the linear resistance vanishes, is impossible in 
two dimensions Iq-fl]. Although the random potential 
associated with disorder can quench the motion of the 
vortex medium as a whole, in two dimensions (in con- 
trast to three) the flow of magnetic flux can be mediated 
by the motion of thermally created point-like defects. To 
the best of our knowledge, only limited experimental ev- 
idence for the absence of the vortex-glass transition in 
two dimensions has been found so far |^,0]. 

Relying on sheet impedance measurements, in this 
Letter we demonstrate that, at typical laboratory time 
scales, the linear response of the two-dimensional (2D) 
vortex medium in ultrathin YBa2Cu307 (YBCO) films to 
a small driving ac field drastically changes with decreas- 
ing temperature. At high temperatures the dissipative 
component of the impedance is frequency independent 
and exhibits a strong (exponential) temperature depen- 
dence which can be attributed to the thermally activated 
motion of single vortices. However, at lower tempera- 
tures, this contribution becomes too small to be observed. 
In this regime the dynamic response is governed by a dif- 
ferent process whose frequency (uj) and temperature (T) 
dependences are found to be consistent with the idea of 
a collectively pinned vortex medium, in which large por- 
tions of the vortex manifold are fluctuating between pairs 
of metastable states in the random potential landscape 
provided by the film microstructure |p|-p^ . 

The crossover to the collective pinning regime has been 
observed over the whole frequency range explored in our 
experiments. Quite remarkably, it is located at a temper- 
ature T*(lo) which turns out to be much higher than the 
estimated melting temperature of a 2D vortex crystal in 
the absence of disorder. Thus, the response we observe 
well below T*{oj) does not originate from a collectively 
pinned vortex crystal, but from a collectively pinned (dy- 



namically) frozen vortex liquid (or, in other terms, from 
a strongly disordered vortex solid). On the other hand, 
a study of the liquid-like response well above T*(cj) as a 
function of the applied magnetic field B reveals that in 
this regime the vortex activation process is controlled by 
a surface barrier mechanism |11| . 

It is worth mentioning that the classical theory of flux 
creep in the collective pinning regime [pPJq] and its ac 



extensions 
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ignore the periodic nature of the vor- 
tex lattice, the vortex medium being treated as an elastic 
manifold interacting with a random potential. Therefore, 
this approach should provide a reliable description of a 
disordered vortex solid (or frozen vortex liquid) as long as 
the frequency of observation is high enough for the mo- 
tion of defects with respect to the driven vortex medium 
to appear to be frozen. 

A possible criterion to distinguish between a genuine 
continuous phase transition to a vortex glass and a sim- 
ple dynamic crossover to a collective pinning regime could 
rely on studies of T*(uj) extending at ultralow frequen- 
cies. If a true phase transition occurs, T*{uj) should satu- 
rate with decreasing u at the (nonvanishing) vortex-glass 
transition temperature, whereas in the opposite case the 
decrease of T*(uj) should be unlimited. In our experi- 
ments T*{uj) decreases by only 5% over three decades in 
frequency, thereby preventing us from drawing any con- 
clusion with regard to the existence of a genuine phase 
transition. Qualitatively, in the frequency range explored 
in our experiments the ac response of ultrathin YBCO 
layers is identical to that of thick films, for which signa- 
tures of a true 3D vortex-glass phase were reported by 
several authors. Thus, from a practical point of view the 
glass- like features we observe below T*{uj) turn out to be 
undistinguishable from those of a genuine vortex glass. 

The experiments were performed on two films YBCO-2 
and YBCO-4, respectively 2 and 4 unit cells thick {i.e, 
with thicknesses d = 2Anm and d — A.Snm), grown 
epitaxially (c-axis oriented) by laser ablation onto (100) 
SrTiOs substrates. Both films were sandwiched between 
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FIG. 1. Arrhenius plots of (a) 1/Lg and (b) R, of YBCO-4 
a,t B — O.OIT for 4 representative frequencies. The dashed 
curve in (a) is the inverse kinetic inductance in zero field. T* 
indicates the crossover temperature of R^ at 99.8kHz. 



nonsuperconducting buffer and cover PrBa2Cu307 lay- 
ers. For comparison, a thick {d — llOnm) YBCO film 
was also studied. Their complex sheet impedance Z was 
extracted from the mutual-inductance change of a drive- 
receive coil system [|2| operated with a conventional lock- 
in detector allowing inductances to be measured with 
a sensitivity of ~ IQpH between 30Hz and lOOkHz. 
The bulk in-plane magnetic penetration depth Xab{T) 
was inferred from measurements of the inverse sheet ki- 
netic inductance I/Lq = uolva[l/Z[B = 0)] = d/^o^afc 
in zero magnetic field and found to fit well, over a 
wide temperature range, a parabolic dependence |1^] 
Kb(T) = A-jf(0)[l-(T/re)2] with Km = 550nm 
and Tc = 48.3X for YBCO-2 and Aab(O) = 340nm and 
Tc = li.lK for YBCO-4. 

The inverse sheet inductance 1/ig = cjIm(l/Z) (which 
measures the degree of superconducting phase coherence 
in the system) and the resistive component of the sheet 
impedance Rz = Re(Z) (which measures dissipation) of 
YBCO-4 in a perpendicular field of O.OIT are plotted 
logarithmically in Fig. 1 as a function of 1/T for a set 
of 4 representative frequencies. As highlighted by the 
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FIG. 2. (a) Log-log and (6) lin-log plots of the 
zero-temperature activation energy as a function of the mag- 
netic field. The dashed lines are fits through the data. The 
dotted line in (a) is a B~^'^ power law. B* indicates the 
2D-3D crossover of the vortex medium in YBCO-4. 



dotted straight lines, at high temperatures both 1/ Lg{T) 
and Rz{T) exhibit an Arrhenius- like behavior, thereby 
pointing to a thermally activated process involving sin- 
gle vortices. In this regime Rz{T) is almost frequency 
independent, whereas l/Lg{T), as shown in the insert 
of Fig. 1, is proportional to uj"^ . These features are in 
excellent agreement with what one expects from barrier 
limited diffusion of non interacting particles. 

The linearity of the Arrhenius plots in the ther- 
mally activated regime implies that the activation en- 
ergy U{T,B) is either constant or linearly dependent 
on temperature. We expect the latter to be the case, 
since U{T, B) should be proportional to the basic en- 
ergy scale of vortex matter, e{T) — 0Q/47r^oA^(,(r)(0o is 
the superconducting fiux quantum and /io the induction 
constant), which varies as {T^. — T) in the temperature 
range of interest just below T^. Then, noticing that the 
temperature dependence of A~j discussed above leads to 
U{T,B) w 2U{0,B)[l - (T/T^)] near T^, our data al- 
low to explore the dimensionality of the vortex medium 
by studying the field dependence of the zero-temperature 
activation energy U{0,B). 

The results for the three YBCO samples studied in this 
work are shown in Fig. 2. While the energy barrier for 
the thickest (reference) film [Fig. 2(a)] obeys a power law 
C/(0, B) ex J5-" with an exponent a « 0.40 ± 0.05 fairly 
close to the prediction a = 1/2 for plastic vortex motion 
in a 3D vortex liquid 0, for the thinnest YBCO-2 layer 
[Fig. 2(b)] U{Q,B) exhibits, over the entire field range 
covered by our experiments, a logarithmic field depen- 
dence signaling 2D behavior. This interpretation is cor- 
roborated by a quantitative comparison with the theoret- 
ical predictions for fiux flow in two dimensions controlled 
by surface barriers |ll[] . In this regime the prelogarithmic 
factor of U{0,B) is of the form Ce{0)d, with C = 1/2. 



Inspection of the slope of U{0,B) for YBCO-2 in Fig. 
2(b) gives dU{0,B)/d\ii{l/B) « 71K, corresponding to 
C ~ 0.45, a value in excellent agreement with that cal- 
culated for the surface barrier mechanism. Note that the 
response of samples of macroscopic size as those studied 
in this work would be dominated by edge barriers only at 
extremely low frequencies, far below those accessible to 
our experiment. However, as evidenced by AFM images 
of the film microstructure, a consequence of the "unit cell 
by unit cell" growth of the YBCO layers is the formation 
of linear defects (steps), related to thickness variations in 
multiples of a unit cell, separating flat islands ~ 0.2^m 
in size. It can be then expected that such steps play a 
role similar to that of sample edges, as they would also 
provide barriers against vortex motion exhibiting, in two 
dimensions, a logarithmic field dependence. 

Additional evidence for the 2D-3D dimensional 
crossover of the vortex medium is provided by a study 
of the activation energy for YBCO-4. As shown in Fig. 
2(a), for this film U{Q,B) exhibits algebraic behavior 
with a ~ 0.58 ± 0.06 at high fields, but crosses over, 
at B* w IT, to a low- field regime characterized by a 
much smaller exponent (a « 0.1) pointing to a logarith- 
mic field dependence, which is indeed demonstrated by 
the lin-log plot of Fig. 2(b). By comparing the activation 
energies for 3D |5[[ll|| and 2D ^^ vortex liquids, it is pos- 
sible to estimate the crossover field as B* « fc(/)o(7/d)^, 
where 7 is the anisotropy ratio and k a numerical con- 
stant of order unity. Using 7 « 1/7 for YBCO, one ob- 
tains B* « IT for YBCO-4 by choosing k « 1/2. 

As shown by Fig. 1 , with decreasing temperature both 
l/Lg{T) and Rz{T) cross over to a regime where their 
temperature dependence becomes much weaker than in 
the activated regime. While the change in behavior in 
l/Lg{T) would be present even if Rz{T) would con- 
tinue to decrease exponentially with 1/T, the crossover in 
Rz{T) points to the onset of a different regime. In order 
to elucidate its physical nature, it is convenient to define 
a crossover temperature T*(uj) and to compare it with the 
melting temperature of the vortex medium. We ad hoc 
identify T*{lu) as the temperature corresponding to the 
maximum curvature of the i?^ (T)-curves in Fig. 1, the 
particular choice of the criterion being irrelevant for our 
conclusions. At i? = O.OIT, T*{llj) varies approximately 
from 67K at the upper limit to 64:K at the lower end 
of the frequency spectrum explored in our measurements 
(the maximum curvature in the l/Lg(T)-curves occurs at 
about the same temperatures). Quite remarkably, these 
crossover temperatures are much higher than the tem- 
perature, Tm = (l)ld/32Tr^V3kBfioXli,iTM) ~ 15K {ke is 
the Boltzmann constant), at which the 2D vortex crystal 
in YBCO-4 would melt due to dislocation unbinding flj] 
in the absence of disorder (notice that the data of Fig. 1 
were taken in the 2D regime well below B*). Thus, the 
picture emerging from the impedance measurements is 
that of a vortex liquid which, well below T*(ci;), appears 



to be frozen at the time scales, 1/w, of our measurements 
and whose response, as is shown below, can be described 
by ac extensions of the theory of collective pinning. 

The low-frequency linear dynamic response of a col- 
lectively pinned elastic vortex manifold has been inves- 
tigated in Refs. P,P| and recently, in a more system- 
atic way, in Rcf. ||10(|. In this approach large (in gen- 
eral anisotropic) portions of the vortex medium (vortex 
bundles) are assumed to be fluctuating between pairs of 
metastable states in the uncorrelated random potential. 
Treating these pairs of states as current-driven two-level 
systems with a size distribution v{r), the vortex contri- 
bution /2i,(T, w) to the specific inductance of the system 
can be shown |^,Q to be, quite generally, of the form: 



lz.{T,Lj)^B' 



drv{r 



2, ,2 






(1) 



where u is the average vortex displacement inside a 
bundle in the direction of the current-induced force, 
r the size of a bundle in the same direction, V its 
volume, E the energy scale which can be associated 
with it, Tc the collective pinning length |15| and r^ ^ 
Tc[{kBT /Uq)\ti{1/u!t)Y/^ a frequency-dependent length 
scale setting the maximum size of the vortex bundles 
which appreciably contribute to the response (C/q and 
T are, respectively, a characteristic energy and a relax- 
ation time related to disorder, x the energy barrier ex- 
ponent). It turns out to be convenient to estimate E 
by considering the compressive contribution E^ to the 
energy. In doing this, one has to take into account 
the dispersive nature of the compression modulus 011(5), 
which in thin films (rf ^ \ab) is always nonlocal. Us- 
ing cii(g) « {B^d/fioXab)l^^ ^^ ^^^ regime of inter- 
est, one obtains [|[ i?c ~ {B'^d/^oXlf^)S'^{u/r)'^, where 
S = V/d ex rr±. Then, setting v(r) oc 1/r^ as imposed 
by the presence of a hierarchical distribution of quasi- 
isotropic {r± ex r) two- level systems iQ, from Eq. (1) 
we find: 

LziT, Lo) « LoiT){l + C' \n[{kBT/U^) ln(l/c^T)]}, (2) 

where C is a numerical constant. Since Lz{uj) depends 
only logarithmically on w, a simplified Kramers-Kronig 
relation flql can be used which leads to: 



Rz{T,iu) « C {7T/2)LjLo{T)/ln{l/LOT). 



(3) 



Notice that in a film with a high density of linear defects 
(steps) acting as a network of strong pinning lines, the 
regime of collective pinning can be realized by vortices 
moving along these linear defects. 

The main features of the response of our very thin 
YBCO layers in the frozen vortex-liquid regime below 
T*{ijj) are well described by these expressions (note that, 
below T*[ijj), Rz <C ojLz and, therefore, Lg « L^). Fo- 
cusing first on the temperature dependence, one sees that 
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FIG. 3. Sheet resistance vs frequency isotherms taken 
at CS/C-intervals between 65K and 69K for YBCO-4 at 
B = O.Oir on a log-log plot. The dashed line is a power-law 
fit of the GS-ft'-data with an exponent « 0.7. The dotted line 
shows the (almost) linear frequency dependence predicted by 
Eq. (3). 



which dissipation could be studied with sufficient accu- 
racy. Considering the fact that this isotherm reflects, at 
a time scale ~ 1/w, the response of a 2D vortex medium 
"on the verge of freezing" rather than that of a "deep- 
frozen" liquid, we interpret the general behavior emerg- 
ing from Fig. 3 and, in particular, the value of the ex- 
ponent extracted from the "coldest" isotherm, as an in- 
dication that Rz{i^) will likely tend to the almost linear 
frequency dependence predicted by Eq. (3) well below 
T*{oj). This should be compared with the response of a 
non-frozen pinned vortex liquid for which one would ex- 
pect a crossover to an w^-dependence at sufficiently low 
temperatures. 
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stages of this work, A. Daridon for the AFM images and 
A. Riifenacht for assistance in the analyis of the data. We 
are also grateful to G. Blatter and V.B. Geshkenbein for 
interesting discussions and useful comments. This work 
was supported by the Swiss National Science Foundation. 



it should be dominated by Lo{T), the logarithmic correc- 
tions entering Eqs. (2) and (3) varying too slowly to be 
of any relevance in the limited temperature interval of 
our experiments. The dashed curve in Fig. 1(a), which 
is simply 1/Lq(T) as inferred from zero-field measure- 
ments, fits nicely the low-temperature (almost frequency- 
independent) l/ig(r)-curves. On the other hand, the 
temperature dependence and the order of magnitude of 
Rz below T*{u!) are compatible with Eq. (3) for any 
reasonable estimate of r (10^^^ — lO^^s). We have also 
found that the response in the frozen vortex-liquid regime 
is only weakly dependent on B in the field range covered 
by our experiments (up to 3r), in agreement with the 
theoretical prediction (derived for B ^ i?c2)- 

Further evidence for the interpretation of the response 
in terms of a collectively pinned vortex manifold emerges 
from the analysis of the frequency dependence of our 
data. As shown in Fig. 1(a), at low temperatures, 
well below T*{uj), l/Lg is almost frequency indepen- 
dent, a behavior consistent with Eq. (2), where the ex- 
tremely slow- varying (double-logarithmic) vortex contri- 
bution can be hardly expected to be noticeable against 
the superfiuidbackground 1/Lq. To discuss the dissipa- 
tive component, in Fig. 3 we show a family of Rz vs uj 
isotherms in a log- log plot. With decreasing temperature 
the isotherms progressively evolve from the frequency- 
independent regime characteristic of the vortex liquid at 
high temperatures to an almost algebraic behavior with 
an exponent ^ 0.7 at the lowest temperature (65 K) at 
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